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Abstract 

Dirac's equation of the electron will be discussed by using quaternions 
as the basis of a new formalism which seems to be very well adapted to the 
problem. The transformation properties of the equations as well as the in- 
variant and covariant [bilinear] constructions of Dirac's theory are developed 
uniformly and systematically. A method of obtaining a covariant formulation 
of the equations using customary tensor calculus also offers itself unequivo- 
cally if we duplicate the Dirac equations. The resulting system consists of 
two parts. The first part corresponds to the coupling of an antisymmetric 
tensor with a vector of the same structure as in the relationship between 
electromagnetic field strength and current vector in Maxwell's equations. 
The other part yields a "feedback" in the form of a reaction of the current 
vector upon the field strength of the same structure as in the case of the 
well-known connection between the vector potential and the field strength 
for the electromagnetic field. 



1 Introduction 

Dirac's equation is based on two viewpoints. On the one hand, it should be a first- 
order linear differential equation and, on the other hand, the iterated application of 
the equation (for the field-free case) should yield the Schrodinger wave equation. 1 

*Editorial note: Published in Zeits. f. Phys. 57 (1929) 447-473, reprinted and translated in [11. 
This is Nb. 1 in a series of four papers on relativistic quantum mechanics [l] |2] |3j |4] which 
are extensively discussed in a commentary by Andre Gsponer and Jean-Pierre Hurni |5 1. Initial 
translation by Josef Illy and Judith Konstag Maskd. Final ttanslation and editorial notes by Andre 
Gsponer. 

^Editorial Note: By Schrodinger wave equation Lanczos means the relativistic Schrodinger 
wave equation which today is usually referred to as the Klein-Gordon equation. 
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The latter aspect in practice implies the invariance under Lorentz transformations 
without stating that the transformations involved should necessarily have a com- 
mon vector analytical meaning. The matrix calculus and the operator method 
have the advantage that they are able to provide a solution to the problem without 
being concerned with the requirements of normal tensor analysis. This advantage 
is diminished by the uncomfortable fact that in a physical theory one has to work 
with quantities which do not allow interpretation by the normal concepts of tensor 
analysis — though these concepts have otherwise proved extremely useful in de- 
scribing natural phenomena. On these grounds, several trials were made to avoid 
the operator method and to bring the Dirac equation into a form which could be 
interpreted in terms of normal vector analytical concepts. C. G. Darwin 2 pointed 
out an analogy with Maxwell's equations in the special case when the mass term of 
Dirac's equation vanishes. On the other hand, Madelung 3 worked out a system of 
equations which can be considered to be a generalization of Maxwell's equations 
and which enables the Dirac equation to be obtained in a natural interpretation. 
For this system, the invariance under Lorentz transformations remains a problem. 4 

In the present paper, the author discusses the problem from a somewhat differ- 
ent point of view. Without heuristically trying to find analogies with the classical 
field equations, he studies the transformation properties of the Dirac equation in 
general by using a formalism which seems to be well-adapted to the problem and 
is based upon the "quaternions" introduced by Hamilton. The quaternion calculus 
has never really been adopted in physics. Uniting vector and scalar multiplication 
into one operation proved to be insufficiently flexible for the purposes of vector cal- 
culation. And later, facing the magnificent unifying and generalizing viewpoints of 
tensor analysis, the operation with quantities based on special properties of three- 
or four-dimensional space had to be withdrawn in favor of a purely component- 
based representation. Nevertheless, there remains the fact that the quaternions 
can offer considerable practical ease in treating the general transformations of the 
Lorentz group, partly because they enable any given Lorentz transformation to be 
described fairly simply, and partly because those constructions, which in terms of 
tensor calculus are to be interpreted as invariant or covariant quantities, are also 
characterized by special and easily describable properties in quaternion calculus. 

The author arrived at the present research through the fact that in his doctoral 
dissertation written 10 years ago 5 he had dealt with the demonstration of the 

2 Proc. Roy. Soc. 120, 621, 1928; Nature 123, 203, 1929. 
3 Zeits. f. Phys. 54, 303, 1929. 

^Remark during revision: I was informed from a kind letter by Prof. Madelung that the 
covariance does exist all the same. The given transformation is linear, but it is not of the usual 
vector analytical nature and thus his system does not come under the viewpoint discussed here. 

5 "The Function Theoretical Relationships of the Maxwell Aether Equations," Publishing House 
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connection between quaternions and Lorentz transformations, and he had pointed 
out how these can be used to formulate the laws of the special relativity theory 
formally in a very simple way, especially when applied to the electromagnetic 
field. Later, he no longer followed this line; however, an accidental observation 
led him to extend these investigations. Initially they were performed with respect to 
Maxwell's equations — later to the Dirac equation by using the same mathematical 
tools, and indeed this course has proven to be correct and natural. The observation 
was that the formerly discussed equation system, which had been considered as a 
generalization of Maxwell's equations, is, interestingly, completely equivalent to 
Dirac's equation for the case where the mass term disappears. The inclusion of 
this term does not present any problem. Not only does it permit a fairly simple 
overview of the transformation properties of the functions, but it also offers new 
prospects concerning the tensor analytic meaning of the Dirac equation and may 
possibly lead to completely new viewpoints. Indeed, the following developments 
are of such almost trivial simplicity, and the direction of the progression is traced 
out so clearly, that one cannot help the feeling of having found a "via regia" to 
penetrate into the essence of Dirac's equation in a way that is more suitable than 
that allowed by the general operator method which is possibly unsuitable for that 
problem. 

Before beginning to discuss our actual problem, we will briefly summarize the 
foundations of the less well-known quaternion calculus and give an outline of the 
most important methods and results of the above-mentioned doctoral dissertation 
(Sections 2, 3, 4). 

2 The quaternions 

By a quaternion Q we understand the combination of four quantities, "compo- 
nents," in the form: 

Q = Xj x + Yj y + Zj z + T Jl . (1) 

The quantities (j x , j y , j z , ji) are four "unit vectors." For the fourth component 6 
we chose the notation / instead of t because we want to reserve t for the real time, 

of Josef Nemeth, Budapest 1919. (Due to the difficulties of the post-war period, the essay was 
published in 50 lithographed copies only.) Inaugural dissertation at the University of Szeged, 
Hungary. Editorial note: This dissertation is now available in paper 1 6 1 and electronic 1 7 1 forms, 
together with two commentaries |8]|9)- 

6 In the pre-relativistic period of Hamilton the fourth quaternion unit had, of course, no rela- 
tionship to physical time and was added as a common scalar to the three spatial units. With regard 
to the applications we are interested in, we have immediately introduced the fourth dimension as 
a time dimension. 
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whereas: 



/ = ict. 



(2) 



Along with the self-understood operation of addition, multiplication is the 
fundamental operation. This is defined by the multiplication of the unit vectors. 
We fix the relations: 

3x3y 3zi 3y3x 3x3yi 3x3 1 3l3x 3x1 

Jx = 3 y — 3 Z = ~3i — —Jl 

The remaining equations are obtained by cyclically interchanging x, y, z. The 
fourth unit vector ji, in the direction of the imaginary "time axis," behaves like the 
ordinary unit. Thus one can also write ji — 1; in other words, j t can be ignored as 
a factor. 

The multiplication is associative but not commutative. Instead of the simple 
commutative law, here we have the law: 



GF = FG (4) 




or 

GF = FG, (4a) 

where the bar means the following: one goes over to the "conjugate" of the 
quaternion, that is to say one gives the spatial components — the space part as we 
shall call them — opposite sign: 

F = -Xj x -Yj v -Zj g + T. (5) 

For an arbitrary number of factors, the conjugate of the product is obtained by 
writing the sequence of factors in reverse order and taking the conjugate of each 
factor. 

It is easy to see that the quantity FF is simply a number (the spatial compo- 
nents = 0). 



3 Four-dimensional rotations and quaternions 

Multiplication of quaternion F by quaternion p in the sense of 

F' = pF (6) 
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can be conceived as a transformation of the line segment F if the quaternion 
is represented as a vector in the four-dimensional space. The matrix of this 
transformation reads: 



/ Pi ~P3 P2 Pl\ 

P3 Pi ~Pl P2 

~P2 Pi Pa P3 

\~Pl ~P2 ~P3 Pi J 



(7) 



(Here the components of the quaternion p are denoted by numbers instead of 
letters.) If we require that the norm of p have the form 



pp = p\+pI + pI + pI = i, 



(8) 



then it is seen that the transformation is an orthogonal one. We shall denote 
an orthogonal transformation of this kind briefly as a "^-transformation." From 
the associative law of multiplication, it follows that the ^-transformations form a 
"group" which is contained in the general orthogonal transformations. 

A second group will be obtained by post-multiplying with the quaternion 
(instead of pre-multiplying): 

F' = Fq. (9) 
This "g-transformation" 7 has the following matrix: 



/ 94 93 -<?2 9i\ 

-93 94 9i 92 

92 -9i 94 93 

V~9i -92 -9s 94/ 



(10) 



The two groups mutually complete each other by supplying in their composition the 
most general group of orthogonal transformations. Thus an arbitrary orthogonal 
transformation can be written in the form 



F' = pFq, 



(ID 



given that: 



pp = 1, qq = 1. 



(12) 



An arbitrary four-dimensional rotation is characterized by six parameters. In 
fact we have six quantities in the two quaternions since their lengths are normalized 
to 1. 

7 In view of the completely different subject matter, this notation should not be confused with 
the customary one used for the generalized coordinates and momenta. 
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A subsystem of the general transformations is formed by the purely spatial 
rotations which do not change the time axis. Then obviously the same relationship 
holds between F and F as holds between F' and F. Since a reversal of the factors 
in (fTTTl yields: 

F=gFp, (11a) 

we must have: 

p = q, q = p. (13) 

Accordingly, purely spatial rotations are characterized by the quaternion at the 
back being equal to the conjugate of the quaternion at the front. 

In the four-dimensional space of reality, one of the dimensions is imaginary. 
Accordingly, an orthogonal transformation cannot contain real coefficients exclu- 
sively, and we should consider both characteristic quaternions to be complex quan- 
tities. On the other hand, they cannot be arbitrary complex quantities. Namely, 
only real Lorentz transformations can be allowed because only these transform 
real (x, y, z, t) into real (x', y', z', t') again. This constraint can be characterized 
as follows: Let us consider the "position vector:" R = (x, y, z, I). Let us change 
over to the conjugate complex quantity which we will denote by an "asterisk" (*). 
As I is imaginary, only the time component of R will change its sign. Hence we 
have 

R* = -R. (14) 
This peculiarity must remain valid in the primed system as well. Hence, if 

R' = pRq, (15) 

then: 

p*R*q* = -qRp (16) 

and so: 

P* = Q, Q* = V- (17) 

Accordingly, real Lorentz transformations are characterized by the quaternion at 
the back being equal to the complex conjugate and "overbarred" value of the 
quaternion at the front. This means that it is enough to give p, and this value 
already determines q. So a real Lorentz transformation can be written in the form 

F'=pFp*. (18) 

In the case of spatial rotations the quaternion at the back must be equal to p and 
therefore it is seen that spatial rotations always belong to real p. 
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4 The Hamiltonian operator 



Let us introduce the following differential operator ("gradient"): 







o 



d d 



3x ~r Jy Q + Jz ~l~ o; 

ax ay az a/ 



(19) 



which we will denote as "Hamiltonian operator." 8 Under coordinate transforma- 
tions it behaves like a vector. 



We shall apply this operation to a quaternion F in the following form: 
VF = 





aZ3x ~ —Jy - -q- z Jz + ^ ) {Xj x + Yj y + Z Jz + T). 



Ox 



Oy 



01 



Let us put: 



VF = 0, 



(20) 



(21) 



so we carry out to some extent the same process in four dimensions, as in the 
complex function theory in two dimensions, if we write down the Cauchy-Riemann 
differential equations by putting: 



0_ 

Ox 



.0_ 

'Oy 



(u + iv) 



Ou 
Ox 



Ov 
Oy 



( Ou 0v\ 
\0y Ox J 



0. 



It was this formal analogy which motivated me at the time to investigate equation 

<eh>. 



If we decompose the above into components, we have the following system: 

OX OT OY OZ _ Q 
Oy 

=0. 



01 
OY 

~bT 
oz_ 
~oi 

ox 



OT OY 

Ox Oz 



OT 
Oy 
OT 



0Z_ 

Ox 
OX 



Oz ~^ Oy 



OX 
Oz 
OY 
Ox 



OY OZ OT 
Ox Oy Oz 01 



=0. 



=0. 



(21a) 



In Hamilton's works we can find only the spatial part of this operator. Strictly taken we should, 
therefore, speak of an "extended Hamiltonian operator." For the sake of brevity, however, we forgo 
this, especially since this extension is quite obvious for the relativistic application of quaternions. 
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These equations are closely connected with Maxwell's equations for empty 
space. There we have: 



(22) 



Subtracting the second equation multiplied by i from the first, we obtain the 
complex equation: 

1 d(H + i£) 
c idt 
in which only the combination: 



ld£ 


-rot ft 




c~dl 


= 0/ 


IdH 


-rot£ 






= 0, 


c~dt 








div£ 


= 0, 




divft 


= 0., 



rot(W + iS) = 0, (23) 



T = H + iS (24) 
occurs. In the same way, we can combine the last two equations in the form: 

div(ft + iS) = 0. (25) 

If we denote the components of T by X, Y, Z (which now also should be 
considered as complex numbers) and write down the equations component by 
component, it becomes clear that the system obtained is identical with equation 
(121 all if we put T = in the latter. 

Therefore equation (|2T1l . which is also Lorentz invariant and also supplies 
the wave equation for all components — just like Maxwell's equations — can 
be conceived as a natural formal extension of the latter. The wave equation is 
obtained by a second application of the V-operation in the following form: 

V(VF) = (VV)F = 0. (26) 

Namely, the V V operator is evidently a scalar and identical with the Laplacian A 
in four dimensions: 

™ a d 2 d 2 d 2 d 2 

VV = A = 1 1 1 . (27) 

dx 2 dy 2 dz 2 dl 2 

5 Dirac's equation for the case of vanishing mass 



If in (12 1 ah we add the second equation multiplied by i to the first, likewise the 
fourth to the third equation, then we have two equations which contain only the 
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combinations X + iY and Z + iT: 




(28) 



Of course, these two equations cannot replace the original four. The eight Maxwell 
equations were halved to four by changing over from real to complex quantities. 
Here, however, X, Y, Z, T themselves are already complex. However, we can 
write two further equations by taking into account the system of equations for the 
conjugate complex quantities and performing the same operations there as well. 
Then obviously we should mark all components with an asterisk and also take into 
account that the term with d/dl changes its sign. Thus for the quantities marked 
with an asterisk, we have the same equations but with an opposite sign in the first 
term. 

Let us introduce the following assignments: 



then it is seen immediately that we have Dirac's equations in the same form as 
they are explicitly written in Weyl's textbook 9 if we neglect the mass term. From 
this we can see that the V operator can substitute fully for the Dirac operator, and 
this has the advantage that the Hamiltonian operator is in closer connection with 
vector analytical quantities. 

We note that the association of Dirac's equations to those of Maxwell has been 
performed by C. G. Darwin 10 in a similar way except that we have to put T = in 
the Maxwellian case. Then the problem arises that ipi and tps are not independent 
of each other; rather, the one is determined by the other. Equation (|2lT) . on the 
other hand, allows exactly as many degrees of freedom as the Dirac equations, and 
they are in fact equivalent to them. 

Before starting to discuss the full Dirac equations, we should like to point out 
an important characteristic of equation (OTT) to which we will return later. The 
equations obviously do not unequivocally prescribe the transformation of quantity 
F under a Lorentz transformation. In fact, without any rotation of the axes, from 

9 Hermann Weyl, Group Theory and Quantum Mechanics (Publishing House of S. Hirzel, 
Leipzig, 1928) p. 171. 
10 loc. cit.; see footnote|2] 



X + iY = 
Z + iT = ip 3 , 



X* + iY* = ^s, 
Z* + iT* = ipi, 



(29) 
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the associative law of multiplication it already follows that a transformation of the 
form: 

F' = Fk (30) 

(with k an arbitrary quaternion) is possible without any change of the system. In 
the case of a Lorentz transformation: 

V = pVp* (31) 
generally one can only say of the transformation of F that: 

F' = pFk (32) 
must hold where k may still be an arbitrary quaternion. 

In the special case of Maxwell's equations, this ambiguity does not appear 
because from the very beginning T = and it has to be demanded that T' = 
holds in the new system as well. Then the rotation imposed on F by transformation 
(l3~2l must also be a purely spatial one and we have: 

k = p. (33) 



We know that the electromagnetic field strength transforms in the four-dimensional 
formalism like an antisymmetric tensor (a six-vector). From this we can draw the 
following conclusion: If we know that a complex quaternion is transformed as: 

pFp, (34) 

then the time component is an invariant of the transformation while the real and 
imaginary parts of the spatial components can be considered as forming a six- 
vector. Of course, we could also give a direct proof for this type of transformation 
of an antisymmetric tensor, but here we shall dispense with this for the sake of 
brevity. 



6 Dirac's equation in the case of imaginary mass 

It is natural to try to introduce the mass term by introducing a term proportional 
to F instead of setting the right-hand side of equation d2TT) equal to 0. However, it 
should be taken into consideration that it is not the repetition of the V-operation but 
the operation V V which gives the wave equation. On the other hand, we can make 
use of the specific characteristic of the four-dimensional space-time continuum, 
that: 

V = -V*. (35) 
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We now make the assumption: 



VF = aiF*, (36) 

where a is any real number. With regard to rule (l35l) it follows from this equation 
that: 

VF* = aiF, (37) 
and so, applying the operation V to the first equation, we have: 

VVF = -c?F. (38) 

Thus we really obtain the required Schrodinger equation, but with the opposite 
sign of the mass term, which would only be possible in the case of imaginary 
mass. Later we shall see that this change of sign is well-founded and cannot 
be eliminated from the equation by introducing numerical factors. Although an 
imaginary mass obviously has no real meaning, it is of heuristic value to discuss 
equation (l36l) briefly. Let us combine the equations in pairs again, exactly as in 
equations (|2"%1) . and let us introduce the quantities tp in the sense of the assignment 
(I2"9"t. Then we obtain exactly the Dirac equation if we take the imaginary value: 

h 

m = ai (39) 

Z7TC 

for the mass m. 

In the case of a Lorentz transformation, the transformation law of quaternion 
F is easily found. Let us set: 

F' = xFy, (40) 

then we have: 

V'F' =p*VpxFy. (41) 

Hence it follows that: 

x*=p*, x = p, y = y* (42) 

must hold if equation (13*61) is to remain valid in the transformed system as well. 
The second condition is equivalent to the first. The third condition states that y 
must be a real quaternion. 

Accordingly, the law of the transformation of F reads: 

F' = pFr, (43) 

with an arbitrary real r, whose length we assume to be normalized to 1. It is 
obvious that no invariant meaning can be assigned to a quantity which behaves 
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in this manner under a transformation. In the case of purely spatial rotations p is 
real; thus we could put r = p and F would be a vector. In the case of a general 
Lorentz transformation, however, such an assignment is impossible. 

This is not remarkable. When writing down equation (BBT) . we make essential 
use of the characteristic of the space-time continuum that one of the dimensions 
is imaginary. In a purely real four-dimensional manifold, an equation of the type 
(l36l) would lose its meaning. However, the concepts of tensor analysis are formed 
in such a way that they never utilize the reality of the quadratic fundamental form, 
and they consider the peculiar values +3,-1 for the index of inertia of actual 
space-time as a coincidence. Here, however, it is just these peculiarities that are 
essentially utilized. 

Even if F itself has no invariant meaning, it can be used to construct quantities 
having such a meaning. First we have an invariant: 

F'F' = pFrrFp = p(FF)p = FF. (44) 

This appears as a complex quantity and is thus equivalent to two real invariants, 
just the two invariants of the Dirac theory. Furthermore the quantity: 

FF* (45) 

forms a vector. Indeed, the transformation law reads: 

F'F'* = pFrr*F*p* = p(FF*)p*, (46) 

since r is real. This vector, the spatial part of which is purely imaginary and the 
time part purely real, is identical with the "current vector of probability" of the 
Dirac theory. 

There exist no further covariant constructions unless a restriction is introduced 
for r, which cannot be made without arbitrariness. 



7 Dirac 's equation in the case of real mass 

It is also possible to arrive at the Schrodinger equation in the following more 
general way. Instead of supposing that VF is proportional to a quantity formed 
from F, we introduce a new quantity G. 

VF = atG. (47) 

If we want: 

VG = 0F, (48) 
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then obviously we have: 

VVF = -a(3F. (49) 

Here the product a/3 can be either positive or negative. By attaching a factor to 
G, the constant on the right-hand side can be made equal or the negative of each 
other, whereby we have reduced the system to two normal forms, namely either: 



VF = aG, 
VC7 = aF, 

or 

VF = aG, 

VG = -aF. 



(50) 



(51) 



Let us first consider the second case. If we change over to the conjugate 
complex quantity in the second equation then we can also replace (BTT) with: 




(51a) 



VF = 

VG* = 

Now we can add these two equations and write: 

V(F + G*) = a(F + G*)* (52) 

or 

V(F + iG*) = ia(F + iG*)*, (53) 

and in this way we have arrived at our former equation combination (l36l) for the 
combination H = F + iG*. 

However, in the first case, corresponding just to a real mass, a unification of 
this type, which should lead to an equation for only one single quantity, is not 
automatically possible. 11 In the following we shall deal with this case. 

It is useful to write our fundamental equation in the form: 

VF = aG* 
VG = -aF*. 



(54) 



Now we can obtain the Dirac equation again if we write the combination (I + 
ill, III + ilV) in both the first and second system. Then we have to make the 
assignment: 

X 1 +tY 1 = fa, X* + iY* = ^ 
Z 1 + iT 1 = tp 3 , Z* + iT* = ^, 



(55) 



"Later we shall see, however, that there exist corresponding combinations here as well if 
quaternions are introduced as factors. 
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and we have for the mass: 



h 



(56) 



m = 



a. 



2tic 



The indices 1 and 2 refer to quaternions F and G, respectively. However, the 
system (154b . consisting of eight equations, is still not exhausted. Now we can 
exchange the roles of F and G and make the assignment: 



A Dirac equation system holds also for the values of a which are independent 
of ip. This equation system is completely identical with the first one, except that 
the signs are opposite in the mass term due to the change of sign in the second 
equation (l54b . These two systems altogether are now equivalent to system d54l . 12 

In the case of a Lorentz transformation d3TT) . the transformation of F and G is 
easily obtained in the following form. It reads: 



where k may be an arbitrary quaternion. 

Here, too, the transformation displays a peculiar uncertainty in ip for which 
no correspondence can be found in the Dirac theory because there the transfor- 
mation of the quantities is unequivocal. This can be explained as follows: For 
a transformation of (X 1 , ...,Ti), (X 2 , ...,T 2 ), after conversion to ip and a, the 
latter quantities will generally be mixed. If, however, we presuppose only one 
Dirac equation, e.g., that containing ip, then only those transformations come into 
consideration in which the quantities a do not appear, i.e., those in which ip are 
transformed into each other. This means that here we distinguish a subsystem 
in equation system d54l) and require that this subsystem should be transformed 
into itself. Thereby the transformation of the quantities ip will be unequivocally 
defined (up to a common factor), but the so-distinguished subsystem will not get 
an invariant meaning, whereas this can be expected of the whole system. In fact, 
this will be confirmed and we shall be able to find a general covariant formulation 
for the extended system. However, before coming to this, it is interesting to have 
a closer look at the Dirac subsystem and its specific transformation characteristics 
— in view of the importance of Dirac 's theory. 

12 A duplication of the Dirac equation can also be found with the Madelung model (loc. cit.) 
where eight complex equations are written down. 



X 2 + iY 2 = a 4 , 
Z 2 + iT 2 = 0-3, 




(57) 




(58) 
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8 Unified derivation of the co variants of Dirac 's the- 



ory 



In accordance with the above, we now presuppose the equation system for ip alone 
and allow only those transformations which convert tp into themselves. We can 
easily show that, under this condition, the indefinite quaternion k of equation (1581) 
can practically be set equal to 1 . We only have to consider that now X + %Y and 
Z + iT may only be transformed in this combination; otherwise, the quantities 
X — iY and Z — iT which can be deduced from the quantities a would also be 
required. It is easy to show that this means the following for the transformation 
matrix: we divide the square of the matrix into four small squares with a horizontal 
and a vertical line across the center. Then in each of these squares the two terms 
along the diagonal must be equal and the other two the negative of each other. 
If we check the matrix © of a ^-transformation for this characteristic, we shall 
see that this condition is actually satisfied. In case of a g-transformation with a 
matrix dTOb . however, this is true for two squares only and it does not hold for the 
remaining two. Therefore, we have to set k\ and k 2 = 1 for the post-multiplying 
/^-transformation and there remains only a multiplication of all with the same 
complex number. 

This transformation is trivial because of the linearity of the equation. It is 
natural to eliminate this by putting k — 1. Then a similarity transformation of 
this kind, which is possible even without any rotation of the coordinate system, 
is excluded and the transformation of F and G becomes unequivocal. However, 
normalization does not go this far in quantum mechanics. There, only the absolute 
value of each complex factor is normalized but not its phase. This is based 
upon the fact that for quantum mechanics it is only the "hermitian operations" 
which play a role — they are not influenced by this phase. Hence, if we put 
k = 1, then we consider a more restricted group of transformations than those 
allowed by quantum mechanics. In this way we surely obtain all co variants that 
have a quantum mechanical meaning and possibly even more. We can proceed 
as follows: We check all covariant constructions of the restricted transformation 
group k — 1 and eliminate afterwards those which are not compatible with the 
mentioned "phase transformation." 

Now let us consider unequivocally the following transformation: 



We do not know the quantities F and G, only those combinations of their com- 
ponents which occur in equations (l55l) . Instead of these quantities we introduce a 




(59) 
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single quaternion H which we assign to the quantities ip in the following way: 

X + iY = V> 4 , X* + iY* = i; 2 A 

z ■ ;r c A . y: ■ rr <.,./ (< 

It is easy to calculate that H is composed of F and G as follows 



Let us put: 
then 



2H = F + G + i(F-G)j z . 

G + F = I, G — F = K, 
2H = I - iKj z . 



(61) 

(62) 
(61a) 



We can easily deduce a differential equation for H . Namely, if we apply the 
operation V to (161 all — keeping in mind that equations d54b hold between / and 
K just as for F and G — then we have: 

2VH = a(K* + il*j z ) = a(-K* Jz + U*)j z 
= ai(I* + iK*j z )j z = 2aiH*j z , 

that is: 

VH = mH* Jz . (63) 

This equation is — as one can easily verify — equivalent 13 to the Dirac equation 
between the quantities ip. u If we had not had in mind an invariant formulation 

13 Of course, either j x or j y could occur instead of j z on the right-hand side, and then we would 
only have to cyclically interchange X, Y and Z during the assignment of the components to ip in 
equations (I60> . 

u Editorial note: Equation d63i . to be called the Dirac-Lanczos equation, is an important result 
of this paper. In Dirac's formulation the 4-complex-component electron field is taken as a 4 x 1 
column vector ^, and the linear operators are 4 x 4 complex matrices. 

In Lanczos's formulation the same 4-complex-component field is a biquaternion H G B = 
Ma(C) = Cl\ : 2 — C^3,o The linear operators are then linear biquaternions functions of biquater- 
nions, which are isomorphic to the algebra of 4 x 4 complex matrices M^C) = Ma(B) = Cl^\. 

In both formulations the operator space has 4 x 4 x 2 = 32 dimensions over the reals. The 
difference is that in the Dirac formulation the field is an abstract 4-component column vector, 
while in the Lanczos formulation the field is directly related to the algebraic structure of spacetime 
because any biquaternion H = s + v is the direct sum of a scalar s and a 3-component vector v. 

Lanczos's formulation is therefore more suitable than Dirac's for studying and demonstrating 
the "classical" aspects of the electron field, and for making comparisons with the Maxwell and 
Proca fields which are usually expressed in terms of scalars and vectors. 

Finally, in terms of Clifford algebras, the Dirac field ^ is a degenerate 8-real-component element 
of the 32-dimensional Clifford algebra Cl±.\ (i.e., an element of an ideal of that algebra) while the 

16 



of the Dirac equation, then we could have started out with this equation from the 
very beginning and it could have formed the basis of our investigation. 

We can also write down for E the transformation: 

H' = pH, (64) 

and it can immediately be seen that: 

H'H' = pEEp = EE (65) 

is an invariant. If we separate the real and imaginary parts: 

EE = A + Bi, (66) 

then we obtain two invariants — these are just the two fundamental invariants of 
the Dirac theory. 

It is also easy to see that the quantity: 

E'E'* = p(EE)p* (67) 

forms a vector. This vector represents the "probability current" in the Dirac theory, 
and its zero divergence was proved by Dirac. 15 

Remarkably, this vector is supplemented by three others which can be obtained 
as follows: Let V be an arbitrary vector, and let us form the following product: 

EVE*. (68) 

We shall prove that this is an invariant. In fact we have: 

E'V'E'* = EppVp*p*E* = EVE*. (69) 

Now let us write this product in the form: 

4 

EVE* = Y^(Hj a H*)V a , (70) 

o=l 

Lanczos field H is any 8-real-component element of the 8 -dimensional Clifford algebra Cl\p. — B, 
which is therefore the smallest algebra in which Dirac's electron theory can be fully expressed. 

For more details see Refs. |5| 1101 ITT1 1121 . where it is shown that the complex conjugation 
operation appearing on the right-hand side of Eq. d63l > is characteristic of the Fermionic character 
of the Dirac field. 

15 Proc. Roy. Soc. 118, 251, 1928. 
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where V a denotes the components of V. However, if ^ B a V a is an invariant, then 
the B a necessarily represents the components of a vector. Consequently we obtain 
a vector with the components: 

B a = H Ja H*. (71) 

In actuality, not only one but four vectors are obtained in this way, since the invari- 
ant (l68l is a quaternion itself and is thus equivalent to four invariants because each 
coefficient of ji remains unchanged in itself. Let us write down the components 
of the four quaternions each in a line, one under the other. In this way we obtain 
a quadratic array containing the components of a vector in each column. 

Now the scheme (I7TT) can be realized in a very simple way. Namely, we can 
consider it as an orthogonal transformation of j a . In fact, it is obviously a real 
Lorentz transformation (since the quaternion at the back of this expression is the 
conjugate complex and barred quantity of the one at the front). Accordingly, 
we form first a p-matrix from H, a g-matrix from H* and multiply these two 
matrices. Applying the resultant matrix to j a , we obtain an array in which only 
the columns and rows are to be interchanged. The components of the individual 
vectors are, however, each contained in columns. Thus the rows of the matrix 
obtained immediately provide the four vectors. The fourth is identical with the 
current vector obtained earlier. 

The transformation matrix obtained is an orthogonal one. Thus we obtain four 
vectors at each point which are perpendicular to each other. One of them is the 
current vector. The length of all these vectors is the same. 16 

The square of the length is the same for all four: 

= (HH)(H*H*) = A 2 + B 2 . 

It is also easy to calculate the divergence of these vectors (which must be an 
invariant). We should then form: 17 

E f£ = £_(a^*) = ( W + * <™*> (72) 
= (Vh)h* - h(Vh)* = u - zr, 

where we put: 

U=(Wh)H*. (73) 

16 This geometrical result is methodologically very interesting, especially in view of the Einstein 
theory of "Distant parallelism" which is just based on such "local n-frames." 

11 Editorial note: In equations M2\ and J73i the parentheses are meant to specify the range of 
action of the differential operators d/dx a and V. 
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Since the differential equation (l63b holds for H, it follows that: 

U = -aij z H*H* = -ai(A - iB)j z , (74) 

and we obtain: 

£ H = - 2aBU (75 > 

This means that only the vector belonging to the third row of the matrix has a 
divergence which differs from zero. 

The remaining three orthogonal vectors — among them the current vector — 
are divergence free. 

We can also write the four vectors obtained in the form: 

S« = HjjT, (76) 

where j, ; is one of the four quaternion units. Only the third and fourth vectors are 
not affected by the phase transformation. The fourth is the previously mentioned 
Dirac current vector. The third, 

= Hj z H*, (77) 

is the one whose divergence differs from zero. 

We can find the tensors of the Dirac theory in a similarly systematic manner. 
Let us form: 

H'jJF = pHj a Hp, (78) 

where j a should be one of the three spatial quaternion units. (The fourth unit 
supplies only the already known invariant.) We already know that a quantity 
which behaves like this under transformation can be considered in its spatial part 
as an antisymmetric tensor if we separate the real and imaginary parts (just as 
we saw in tile case of the electromagnetic field strength). The time part gives an 
invariant. This time part, however, is zero for all three constructions and therefore 
gives nothing. 

The phase transformation is only compatible with the quantity formed with j z : 

Hj z H, (79) 

which is thus the only quantum mechanically allowed one. This antisymmetric 
tensor was introduced by C. G. Darwin. 18 

18 Proc. Roy. Soc. 120, 621, 1928. 
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There is another possible way to obtain a tensor. Using the vectors U and V 
we form the invariant: 

H'U'V'H' = HppUp*p*VppH = HUVH. (80) 

Let us write this in the following form: 

J2(HjJ u H)U,V u , (81) 

then we have a bilinear form the coefficients of which must be the components of 
a tensor: 

T lk = Hjd k H. (82) 

Here, too, we obtain four tensors simultaneously since the components appear as 
quaternions and a tensor component can be separated from each unit vector. 

We could just as well have started from the invariant: 

H*UVH* (83) 
and, correspondingly, we would have arrived at the tensor: 

T: k = H*j t j k H*. (84) 

With this we have listed all those covariant constructions of zero, first and second 
degree for the restricted group k — 1, which are built quadratically from the 
fundamental quantities. 

The covariants of the Dirac theory have already been discussed in the litera- 
ture. 19 However, the uniform development outlined here may outdo the method- 
ological persuasive power of other descriptions by its clarity and simplicity. 



9 Failure of the current vector with respect to strict 
covariance 

Once again surveying our train of thought, we can state the following: We started 
from a larger equation system and found that the transformation of the functions 
is not unequivocally determined by the equations. However, we required that a 
given subsystem of the system transforms into itself and this requirement enabled 
us to eliminate the uncertainty of the transformation. 

19 See especially J. von Neumann, Zeits. f. Phys. 48, 868, 1928. 
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However, if we only wish to investigate the Dirac equation, we could just as 
well have taken equation (l6"3t as a basis from the very beginning. (We have derived 
this equation from the larger system for a certain combination of F and G.) We 
could have performed the assignment (i60l) and could have shown that equation 
(l63t is indeed equivalent to the Dirac system. If we now take this path and want 
to calculate the transformation of the function H immediately from this system, 
then we shall experience a peculiar discrepancy with our last results, and this is 
probably not without significance. Let us write the transformation of H in the 
form: 

H' = pHk, (85) 
then we obtain the following condition for quaternion k: 

Jzk = k* Jz . (86) 

This condition by no means represents a strict limitation on k that would allow, for 
instance, only a trivial similarity transformation. Rather, condition (l86l) formulates 
the requirement that the x and y parts of k must be purely imaginary and its z 
and / parts purely real. Now if we also perform a natural length normalization, 
which excludes just the trivial similarity transformation, then a very important 
3-parameter group of transformations will still remain which certainly has nothing 
to do with the "phase uncertainty" of the ij). Rather, one can observe that in the 
case of this transformation which is possible even without any rotation of the axes, 
the new ip can be expressed not only in terms of the previous ij) but also in terms of 
their conjugates tp*. As tp only plays the role of auxiliary quantities, no objective 
reasons can be found for excluding these transformations, even more so because 
these are absolutely normal transformations for the (X, Y, Z, T) — without the 
appearance of the conjugates. No reason can be found why we should not introduce 
these quantities instead of ip as fundamental quantities, though obviously the 
Hamiltonian operator is more closely connected with tensor analytical quantities 
than is the Dirac operator. If, however, we allow these transformations, then a 
large part of the covariants set up will be lost. Only the invariant HH and a vector 
will remain. However, this is not the current vector with zero divergence but vector 
whose divergence proved to be different from zero. In fact, we obtain: 

H'jJI* = pHk Jz k*H*P* = pHjXF, (87) 
because from (ISBT) it follows that: 

3z = k*j g k = hjX- (88) 

For the current vector, we cannot prove the invariance under this /^-transformation, 
Thus is the only vector which really seems to be strictly covariant. In the 
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case of the other constructions (especially with the current vector), the covariance 
is the result of what is, after all, an arbitrary constraint. 20 



This circumstance seems to suggest that the Dirac equation should be consid- 
ered as a component of a larger system, instead of a system closed in itself. 

10 Covariant formulation of the doubled Dirac equa- 
tion 

Let us now consider system (l54b as a whole and let us study its transformation 
properties. As we have seen, the system is equivalent to two simultaneous Dirac 
equations which contain two independent groups of quantities: tp and a. However, 
from the standpoint of the whole system, such a decomposition would be rather 
constrained and unnatural and of no advantage even as a mathematical aid, since 
it would amount to separating out a subsystem which is not internally preferred. 
We want to refrain from doing this and to be guided rather by the point of view 
that a covariant meaning should be attributed to the quantities F and G contained 
in the equations. This principle will lead us to remove the uncertainty of the 
transformation which lies in the arbitrariness of the quaternion k, whereas until 
now we attained this just by separating out a subsystem of equations which is not 
covariant in itself. 

Let us once again write down the transformation equations of our functions F 
and G: 



where k may be arbitrary and is only subject to the always possible, but natural, 
normalization kk = 1. Before we come to the determination of k on the basis 
of the covariance principle, let us note that certain invariants and covariants are 
possible even with k left arbitrary. 

Namely, first the two invariants: 



Remark during revision: The author had not suspected that the transformation properties could 
change due to the introduction of the vector potential. Actually, this is exactly the case; thus the 
special group of transformations found here is lost with the extension of the system by the external 
field. The objection that was made against the customary transformation theory of the ip functions 
is thus invalid. For more details, see the paper "On the covariant formulation of Dirac's equation" 
to be published in the near future. (Editorial note: See Ref. |2|.) 




(58) 



FF, GG, 



(89) 
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which are equivalent to four real invariants. 



Then the two vectors can be formed: 



FG* 



GF\ 



(90) 



so that, for example, 



F'Gf* = pFkkG*p* = p(FG*)p*. 



(91) 



The two vectors of (|9"0"1) are, however, not independent of each other; rather, the 
second is equal to the conjugate complex and overbarred quantity of the first. 
Consequently, it is enough to consider FG alone. 

We might guess that this vector — i.e., its real and imaginary parts — should 
correspond to the current vector, as it represents the analog construction of HH . 
However, if we calculate its divergence using a method quite similar to that applied 
in (1721 . we shall find this is not confirmed. Instead we shall have: 



This finding could be used as a negative factor against the interpretation sug- 
gested here. In fact, this objection would not be justified. This is because this 
vector, distinguished by the fact that it remains covariant even under an unlimited 
transformation of the equations, does not correspond at all to the current vector of 
the Dirac equation but to vector (177b . the divergence of which does not vanish even 
there. However, the Dirac current vector cannot unequivocally be made to corre- 
spond to a covariant formulation, as it is not characterized by invariant properties, 
and its covariance is due rather to an unnecessary limitation of the transformations 
— as was found in the last section. 

In our formal introduction, we were able to find a tensor analytical interpreta- 
tion for two types of quantities. If k = p*, then we have to deal with a vector. If 
k = p, then we can consider the construction as a non symmetric tensor — with its 
time component an invariant. Strangely enough, whichever choice we make both 
constructions occur. The difference is only that F and G exchange their roles. 
Now we want to choose: 



By this choice we put F as an antisymmetric tensor and Gas a vector. Now only 
covariant quantities occur in our equation. Also, both systems of equations are 
already covariant in themselves, not just the whole system. 

Let us examine now the first system. This system seems to be closely related 
to Maxwell's equations of the electromagnetic field. Here we are not thinking of 



div(FG*) = S[VFG*} = a[(FF) + {GG)*}. 



(92) 



k = p. 



(93) 
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the vacuum equations but of the complete equations with the current vector on 
their right-hand side. Namely, Maxwell's equation containing a current vector can 
be brought to the following form by the complex combination (l24h : 

VF = S*, (94) 

where S means the current vector considered as a quaternion. This is, however, 
just our first equation which is thus given a very simple interpretation with the aid 
of a classical analogy. Now we should only put F 4 = for Maxwell's scheme 
but the occurrence of this term means only that the gradient of a scalar appears on 
the left-hand side. In addition, the complex character of our vector G should be 
interpreted so that not only an electric but also a magnetic current must be taken 
into consideration. 

The emergence of the second equation is essentially new and unknown to the 
classical theory. This equation means that there is a "feedback" between current 
and field intensity, a reaction of the field to the current. The appearance of this 
equation should be looked upon as the actual effect that led to the discovery of the 
Schrodinger equation. 

The second equation also has a simple covariant meaning. An equation of this 
type is likewise well-known to us from electromagnetic field theory. If, namely, we 
write down the electromagnetic field by means of a vector potential $j, regarding 
the $j vector as a quaternion, then we obtain the equation: 

V$ = —F*, (95) 

which is just our second equation where we only have to consider that in our 
case $ is to be treated as a complex quaternion. This implies that in addition to 
the customary "rotation construction," the "dual" construction appears as well, in 
accordance with the occurrence of a "magnetic vector potential." 

Thus we can convert our equation system into the common language of physics 
either with the customary vectorial symbols or by using the terminology of tensor 
analysis. For this we only have to apply the corresponding equations of the elec- 
tromagnetic field and substitute the corresponding quantities. The only difference 
is that the continuity condition for the current cannot be inferred from the general 
structure of the equations because of the occurrence of the redundant invariant 
which is not contained in the theory of Maxwell. This invariant is responsible 
for the fact that the disappearance of the divergence for the current vector is not 
a necessary consequence of the equations and that the equation system does not 
contain any internal correlations. (Absence of identities.) 

For the sake of clarity, we shall write down the resulting equations explicitly, 
in tensor analytical form. Then the following quantities appear: 
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An antisymmetric tensor: F ik j= —F ki ("electromagnetic field intensity"). Let 
the "dual" tensor be denoted by F ik . 21 

Two vectors: Si and M { ("electric and magnetic current"). 

Two invariants: S and M. 

We obtain altogether 16 equations which can be divided into two groups 
(A and B). Group (A) contains two vector equations; group (B) contains an 
antisymmetric tensor equation (six equations) and two scalar equations. 

With the customary symbolism of tensor calculus, the equations read as fol- 
lows: 



(A) 



0$ 
dx k 



■ dS 


^Fi, 
dx^ 


< 

dM 


qK 


k dxi 


dx^ 


dS k 

+ 

dxi 


(dMi 
\dx k 



- aSi, 

-- aMf, 

~dM k 
dxi 



aF ik , 



dS» 
Ox,, 

dx„ 



aS, 



aM. 



(96) 



Here the constant a is defined as 



a 



27rmc 
h ' 



(97) 



The whole system, just like the Dirac equations, displays a strong symmetry. 
However, it is questionable whether all quantities have a real meaning. So we 
might guess that the two scalars S and M do not appear in reality. The first 
group of equations would then be fully equivalent to Maxwell's equations. The 
continuity condition for the current would be a consequence of the equations. 
Despite this specialization, the previously mentioned difficulty experienced with 

21 The dual assignment should be performed according to the following scheme: 

F12 = 1F34, ...,F\4 = iF 2 3 7 ... 
(The dots indicate cyclic interchanging of 1, 2, 3.) 
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the assignment made by Darwin does not appear here. Namely, the two groups 
of -0-quantities (ip and a) are assigned not to the field intensity alone but to the 
combination of field intensity and current vector simultaneously, and both types 
of quantities appear in each group [as can be seen from equations d55t and (l57l) l . 
This is why now there is no algebraic relationship between the quantities ip despite 
the specialization. 

Moreover, a further heuristic aspect can be derived from Maxwell's equations, 
namely, for the "magnetic current vector" Mj which vanishes there. We must, 
however, take into consideration that here we can by no means dispose of quanti- 
ties so freely as in common field theory. Namely, due to the double coupling, there 
are no dependent and independent quantities here and the number of equations 
is just as large as the number of unknown quantities. It is just the characteris- 
tic feature of quantum mechanics that it operates with homogeneous equations 
("eigenfunctions"), and the extraneous functions (energy, vector potential) appear 
not as "right-hand sides" of equations but as factors. If we put in our system 
certain quantities equal to zero, this would lead to overdetermination: the number 
of equations would be larger than that of unknown quantities. Then, accordingly, 
such a zero assumption should be compensated for by the omission of equations to 
avoid possible contradictions — unless there are a number of identities, whereby 
the redundant equations would seem to be pure consequences of the others. 

From a mathematics and aesthetics viewpoint, this redundancy would be of 
considerable benefit. It cannot be denied that the strong symmetry of Maxwell's 
equations with respect to duality between electric and magnetic field intensities 
exists only as long as the equations are written down with the symbols of three- 
dimensional vector analysis. With respect to the four-dimensional tensor analytical 
description, there is an internal difference between the two systems of Maxwell's 
equations. In reality, this difference manifests itself in the absence of magnetic 
current. In one of the systems, a normal divergence appears, whereas the other 
system contains the divergence of the "dual" field intensity. The tensor analytical 
formulation is faultless, yet it makes use of an accident in the four-dimensional 
space: the fact that in the case of a second-rank tensor, the "dual" tensor will be 
exactly of second-rank again. In principle, the emergence of this circumstance is 
logically not really satisfying. 

Now, if we put Mi = for the magnetic current in our equations (I96T) . then 
we must simultaneously omit only the equations containing the divergence of 
the dual tensor from group (A) which corresponds to Maxwell's system. These 
equations are then consequences of group (B) and no longer belong to the system 
as determinative elements. Then the "dual" construction would no longer occur in 

22 loc. cit.; see footnote|2] 
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the system. We would have ten equations for ten quantities, on the one hand the 
four Maxwell equations: 

BF 

~rr^~ = aS t , (98A) 



and, on the other, the six feedback equations: 

dS 1 _dSk _ 
dxk dxi 



aF ik . (98B) 



The divergence equation for charge is now a consequence of the equations as are 
the missing four Maxwellian equations. 23 

However, this system has become so non-symmetric that its immediate con- 
nection with the Dirac equations seems to be questionable, though it is also of 
first-order and the Schrodinger wave equation is obtained for each component 
here as well. Indeed, we arrive at the Dirac equation again if we complete our 
system by the second-order Maxwell equations and by the divergence equation 
for current — which really hold as identities even if they do not add anything 
new to the system. The eight Dirac equations for the quantities ip again form a 
subsystem of our whole system and contain somewhat less information than our 
ten equations. However, here we already encounter a constraint in the tp quantities. 
For we can see from our assignment (l55l) that ipi = S z — iS t becomes a purely 
real quantity — since the G vector is now no longer complex but real in its space 
component and imaginary in its time component. The remaining ip quantities do 
not show a similar specialization. 

It is obvious that the question as to whether a real meaning can be attributed 
to the above interpretation cannot be answered until the equations are completed 
for the presence of an external field. Our goal did not take us so far. We raised the 
question as to whether the Dirac theory can be formulated by exclusively using 
field theoretically meaningful quantities. We have found a way which was almost 
unequivocally compelling and indeed led to the expected results. It resulted in an 
internally highly consistent covariant system which permits new prospects in some 
respects. Also, a certain arbitrariness — which, as we have shown, is inherent in 
the Dirac transformation theory and suggests an underlying reason — is eliminated 
here, since here any uncertainty is avoided in the transformations. 

Should this chosen way prove unsuccessful, it would most likely be hopeless 
to expect a field theoretical background for the Dirac theory. Our developments 
reveal little doubt that this background can lie only in the direction found here, if 
it exists at all. 



^Editorial note: Eqs. ( 19 8 Al l and J98BI are the correct wave-equations for a massive spin 1 
particle, to be rediscovered by Proca in 1936. For more details, see section 11 in 1 5 1. 
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In formal respects, one can consider it a gain that a method has been worked out 
which enables the transformation properties of the Dirac quantities to be described 
in a very unified and transparent way. 

Berlin-Nikolassee, July 1929. 
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